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Abstract 

We take a fresh look at high energy radioactive nuclei data reported in the 90's and at the 
positron data recently reported by PAMELA. Our aim is to study the model independent 
implications of these data for the propagation time scales of cosmic rays in the Galaxy. 
Considering radioactive nuclei, using decaying charge to decayed charge ratios - the only 
directly relevant data available at relativistic energies - we show that a rigidity independent 
residence time is consistent with observations. The data for all nuclei can be described by 
fs.i = (tj/100 Myr) 0,7 , where f s ,i is the suppression of the flux due to decay and U is the 
observer frame lifetime for nucleus specie i. Considering positron measurements, we argue 
that the positron flux is consistent with a secondary origin. Comparing the positron data 
with radioactive nuclei at the same energy range, we derive an upper bound on the mean 
electromagnetic energy density traversed by the positrons, Ut < 1-25 eV/cm 3 at a rigidity of 
1Z = 40 GV. Charge ratio measurements within easy reach of the AMS-02 experiment, most 
notably a determination of the Cl/Ar ratio extending up to 1Z ~ 100 GV, will constrain the 
energy dependence of the positron cooling time. Such constraints can be used to distinguish 
between different propagation scenarios, as well as to test the secondary origin hypothesis 
for the positrons in detail. 

1 Introduction 

The sources of primary cosmic rays (CRs), the correct description of CR propagation and the 
mechanisms of CR trapping in and escape from the Galaxy are all essentially unknown [Tl [2] [3] . 
Given the list of open questions, it is remarkable that very few detailed model independent 
analyses exist in the literature. In contrast, models of homogeneous diffusion have become a 
standard, where the general approach is to fit the parameters of the model to B/C and low energy 
10 Be/ 9 Be data. It is difficult to extract generally applicable information from such studies. 

In this work we present a model independent analysis of the propagation time scales of CR 
nuclei and positrons. These time scales are related to the CR residence time in the Galaxy, a 
key unknown of CR propagation. Motivated by the prospects for improved measurements in 
the near future [H El [5J [7J [5] , our analysis demonstrates that it is possible to extract significant 
quantitative information from CR measurements under general assumptions, without committing 
to any particular propagation model. 

The main data we use are based on the decaying charge to decayed charge ratios Be/B, 
Al/Mg and Cl/Ar [3], measured to relativistic energies jTD] and studied more than a decade 
ago by [HI . Earlier analyses of CR propagation time scales in the context of diffusion and 
leaky box models can be found e.g. in [T31 [TH [T^l EH [TTJ HE] and references therein. Being 
model dependent, the values deduced in these studies for the CR residence time differ by order 
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of magnitude, as well as by their physical interpretation. In addition, these studies focused on 
low energy isotopic data, where various theoretical complications arise and where the limited 
dynamical range of the experiments precludes a direct inference of the energy dependence of the 
residence time. Apart from the work of [UJ H2], the high energy charge ratio measurements 
we analyze here were recently considered in [Jj|J [50] in the context of fits to the parameters of 
diffusion models. 

Before proceeding to describe the plan of this paper we first emphasize the key concepts and 
approximations. A central tool in this work is the CR grammage, the mean traversed ISM column 
density experienced by CRs (TJ. At high energy, where propagation energy gains/losses can be 
neglected, the local flux of any stable secondary nucleus can be calculated reliably as [21] [22] 

c ~ 

Js = 7- A csc Qs, (1) 

where A csc is the CR grammage and Qs is the local net production (including spallation losses) 
per unit column density of ISM. Eq. (fTJ is an empirical relation; any acceptable model of CR 
propagation must reproduce it. 

In our study we identify the effects due to decay in radioactive nuclei and due to radiative 
losses in positrons, by using Eq. (fTJ to compute the flux that would be expected had there been 
no decay and loss and relating it to the flux observed in practice. Doing so entails some approx- 
imation. For Eq. (TJJ) to apply, the CR rigidity should remain unchanged during propagation, 
including fragmentation and decay events. In contrast, the decaying charge and decayed charge 
samples, reported at a given kinetic energy per nucleon, exhibit a spread in rigidity. This spread 
is more pronounced for lighter elements, the maximal change occurring in the decay 4°Be — ^ B 
in which the rigidity of the nucleus drops by 20%. Similar effects occur in certain spallation 
reactions and in the production of positrons and antiprotons, which do not inherit the rigidity 
of the primary CR. Due to these effects, typical calculations based on Eq. (fTJ are accurate only 
at the ten percent level or so [33] [23] . 

In addition, various propagation effects cause energy change for CR rigidities 72 < 10 GV. 
Striking examples of such effects are the latitude dependent geomagnetic cutoff of the earth 
and the time dependent (charge dependent and independent) solar modulation [3] [3J. Other 
mechanisms likely operate in interstellar space. Since the theoretical understanding of these 
effects is poor, we limit our analysis to CR rigidities 72 > 10 GV. 

The plan of this paper is as follows. In Section [5] we study the high energy radioactive 
nuclei data analyzed by jll] . We present the data in a new way, which makes the existence of an 
underlying time scale in the problem - the CR residence time - manifest. Our presentation of the 
data demonstrates that the charge ratio analysis of [TJJ does not suffer from gross systematics. 
We show how a combination of data from different nuclei species can be used to measure the 
functional form of the suppression due to decay. We show that a CR residence time which is 
mildly dependent or constant as a function of rigidity is consistent with the data. 

In Section [3J we consider the PAMELA positron measurements [U [35] . We argue that the 
positron flux is consistent with a secondary origin [22] . Assuming secondary origin we compare 
the effects of cooling in positrons to the effects of decay in nuclei, extracting bounds on the cooling 
time of the positrons. These bounds can be used to test in detail the production mechanism of 
the positrons, and have implications for CR propagation time scales in general. 

In Section 0] we summarize and discuss our results. We point out that new measurements 
expected in the near future from the AMS-02 experiment [5] can significantly improve our quanti- 
tative understanding of CR propagation time scales in a model independent manner. Our simple 
analysis can readily be repeated for the new AMS-02 data. 
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In Appendix[X]we give two propagation model examples to illustrate some general arguments 
made in the paper. 



2 Lessons from radioactive nuclei data 

The measured effect of radioactive decay in nuclei receives different interpretations, depending 
on the adopted propagation model. At low energies E < 400 MeV/nuc, measurements of the 
isotopic ratios 10 Be/ 9 Be, 26 A1/ 27 A1, 36 C1/C1 and 54 Mn/Mn yield, within a leaky box model 
interpretation, a residence time t csc ~ 15 Myr |15j . Using the same low energy data, diffusion 
models yield a residence time t CS c ~ 100 Myr 18 . Since the isotopic ratios are only known at 
sub-GeV energies, they give no direct information about the energy dependence of the residence 
time and of the suppression of the flux due to decay. 

In contrast, the decaying charge to decayed charge ratios Be/B, Al/Mg, Cl/Ar and Mn/Fe 
have been measured up to energies E ~ 16 GcV/nuc (rigidities 1Z ~ 40 GV [TO])- Ref. [UJ 
interpreted these charge ratios using three different models: a leaky box model, a diffusion 
model and a diffusion-based Monte-Carlo calculation. From the Monte-Carlo analysis, a CR 
residence time t esc ~20-30 Myr was deduced. Given that the usual leaky box and homogeneous 
diffusion models assume t csc oc A csc and that A csc ~ 1Z~~ 5 , it is interesting to note that |11] 
quotes a single value for the CR residence time. Since the data spans about a decade in rigidity, 
the deduced value of t osc would have changed by a factor of three in this range had the naive 
expectation been applicable. 

The charge ratios Be/B, Al/Mg, Cl/Ar and Mn/Fe are predicted in a leaky box model to 
vary by about 20%, 30%, 50% and 40%, respectively, between no decay and complete decay. 
While this estimate may not apply in general, it makes clear that systematic uncertainties are 
of major concern [T51 [TB] . Measurement errors were discussed in [TT] . Here we emphasize three 
types of systematic uncertainties which are most relevant for our analysis. First, uncertainties 
in the fragmentation cross sections are of order ten percent |26j . Second, a potential source of 
uncertainty is the existence of some primary component in the 26 Al, 36 CI and 54 Mn fluxes (2"Tll28| . 
Third, as noted in the introduction, the use of the CR grammage relation, Eq. (jTJ , introduces 
additional theoretical inaccuracy of order ten percent. All three types of bias are difficult to 
quantify precisely. A dedicated analysis, beyond the scope of this paper, would be required to do 
so. In the rest of this section we show that despite this caveat, the results of [TT] contain clear 
physical information, which can be used to constrain the CR residence time directly, in a model 
independent manner. 

Following [22], we define the suppression factor due to decay of a radioactive nucleus specie 

i as 

f.i = ■ (2) 

' n. y 

The suppression factor / S) j can be derived from the surviving fraction fa reported in |11) . using 
the formula 

fs,i = 7 ryi (3) 

l + ^A csc (l-//) 

where <Ji is the total fragmentation cross section per ISM nucleon. Using the suppression factor 
f 3t i rather than the surviving fraction fa allows us to (i) isolate the effect of decay from that of 
spallation losses and (ii) later on, compare directly the effect of decay of radioactive nuclei to 
the effect of energy losses of positrons. 
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It is instructive to examine the suppression factor due to decay vs. observer frame lifetime. 
This is done in Fig. [TJfor 10 Be, 26 Al and 36 C1, using the fragmentation cross sections and rest 
frame lifetimes from Table [T] and the CR grammage from [35]. We omit 54 Mn, the lifetime 
of which is uncertain [15], and adopt a 15% uncertainty estimate for the CR grammage and 
the fragmentation cross sections. A number next to each data point denotes the rigidity at 
that point in GV. We make the following observations: (i) The suppression factor for all three 
nuclei species exhibits a clear correlation with lifetime, (ii) Where available at the same lifetime, 
/s, 10 Bc ~ /s. 36 ci an d f s , 26 A\ ~ fs. 36 c\ to better than a factor of two. 

Since the rest frame lifetimes of 10 Be and 36 C1 ( 26 A1 and 36 C1) are different by a factor of five 
(three) , data points of equal lifetime are separated by similar factors in rigidity. The approximate 
continuity of across nuclei specie thus hints that the rigidity dependence of the residence time 
is not strong. To assess this point quantitatively, we make some rather general assumptions. 



reaction 


t [Myr] 


a [mb] 


i°Bc -> l°B 
? 3 A\ -+ ?«Mg 
3 fCl -> ?jAr 


2.18(0.09) 
1.31(0.06) 
0.443 (0.003) 


210 
410 
520 



Table 1: Decay lifetimes and fragmentation cross sections for unstable radioactive isotopes. 
Decay lifetimes refer to pure j3 decay channels [TSJ [17] . Fragmentation cross sections refer to 
spallation on H. We adopt a 15% uncertainty estimate for the cross sections and correct them 
for an ISM composition of 90%H+10%He, using the prescription in [3D]. 



It is natural to expect the residence time t csc to be a function of rigidity, 

^csc ^escV^y) (4) 

and f St i to be a function of {U/t esc ) where U is the observer frame lifetime. Limiting the discussion 
to high energy CRs with 72. > 10 GV, the data spans about a decade in lifetime. In this range, 
we assume an approximate power law form 

\ '•CSC / 

Two comments are in order: (i) In general, a depends on (ii/i eS c)0 In writing Eq. (0, we 
assume that this dependence is mild. We will see that this assumption agrees with the data. In 
the future, with data sampled more densely and extending to higher energy, the method we use 
can be straightforwardly extended to account for varying a. (ii) The precise physical meaning of 
t csc appearing in Eq. ([5]) is model dependent. What is relevant for our discussion is that Eq. ([5]) 
is applicable to any nuclei specie (and later on, to positrons) with the same time scale t csc . In 
Appendix lAl we demonstrate how Eq. ([5]) arises in specific models. 

Under assumptions Eqs. (@][5]), comparing the suppression factor for two relativistic nuclei i 
and j at the same rigidity TV we have 

lo HOTor alog l^J- (6) 



In the limit t% 2> tcsc we must have a — > 0. Otherwise when t; < t eS c we must have a > 0. Hence, a decreases 
with increasing (ti/t esc )- 
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lifetime [Myr] 



Figure 1: f s ^ vs. observer frame lifetime. Numbers near data points denote the corresponding 
rigidities in GV. 



Thus we can extract a, regardless of the particular form of t csc . A determination of a from Eq. ([6]) 
requires nuclei with significantly different rest frame lifetimes, since Aa cx 1/ log (Ti/rj). In the 
left panel of Fig. we plot values of a, obtained by applying Eq. ^ to pairs of /io Be , f36 Cl and 
/ 26 Ah/ 36 ci measured at approximately the same rigidity. The result is consistent with constant 
a. Assuming constant a we combine the ten data pairs at 1Z > 10 GV, obtaining 



0.7 ±0.2. 



(7) 



The error estimate in Eq. ([7]) is twice that arising from the errors in f s ^ and in the nuclei lifetimes. 
The inferred value of a is shown as a gray band in the left panel of Fig. [2] 

Using Eq. Q and Fig. Q] we estimate t csc ~ 100 Myr in the range K ~ 10-40 GV. To 
make contact with the correlation depicted in Fig. [TJ comparing the suppression factors for two 
relativistic nuclei i and j at the same observer frame lifetime t' we have 



fs,j (t ) \ ^csc 



where 



eTZi 



t'Aj 

Ti Zi 



Aj ZiTi 
A~Z~V 3 



(8) 



(9) 



Given that Tl3f>c\ ~ ~ 3 7?.26 A1 and the preferred range a > 0.4, we verify that Fig. [TJ 

constrains the rigidity dependence of the residence time. In particular, assuming a power law 
form 

tcsc = t esc< [ — ^ J , (10) 
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Figure 2: Left: Values of a extracted from pairs of / Sj iob C 7 fs, 3B c\ an d /s, 26 Ab /s, 36 ci measured at 
approximately the same rigidity. Right: The residual rigidity dependence of f s j normalized to 
f s ,i = (ti/t esc ) a , with t csc = 100 Myr and a = 0.7. 



we can extract aS from 



log 



fs,i (f) 
fsj (*') 



aS log 



A 3 Zj n 
A t Z 



3 '3 



(11) 

9 and 10 GV 
(12) 



The two pairs of /iob i/ 36 ci an d / 26 ai>/ 36 ci data points at minimal rigidities 
give 

a(5 = 0±0.6 (ftwlOGV), 

were the error estimate is twice that arising from the errors in / Sj j and in the nuclei lifetimes. 

In the right panel of Fig. [5] we plot the residual rigidity dependence of f a + normalized to 
f s l = (t 2 /100Myr) u , inferred from Eq. Q and Fig. [U We find no significant residual bias. 
Thus, Fig. [5] contains the main result of this section: the high energy charge ratio measurements 
are adequately described by Eq. ([5]), implying nontrivial constraints on a and consistent with a 
rigidity independent CR residence time. We make some final comments as follows: 

• As seen from Eq. ©, a can be determined independently of the residence time from the 
normalization of f s .i for different nuclei species. There are five CI data points with 1Z > 10 
GV and about twice that number of Be and Al points with comparable rigidities. Therefore, 
the current data gives a fair determination of a ~ 0.7. 

• In principle, once a is determined, the residence time would also be measured. In practice, 



since for some data point t e 



U/f, 



l/a 



the measurement of t csc is sensitive to the lower 



range of a. Because of the limited range in rigidity, S and i e sc,o are correlated. Thus large 
values of i e sc,o and correspondingly large negative 5 are still allowed in conjunction with 
small values of a. 

• Negative values of 6 < —1 arc in some tension with observations. This is relevant to the 
interpretation of the PAMELA positron data, discussed in the next section. If 6 < —1, 
the positron suppression factor can rise with energy even under the simplest assumption 
of constant radiative losses. 
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• a is significantly larger than zero at least up to t, ~ 30 Myr and TZ ~ 30 GV. Thus, AMS-02 
data [5] (in particular CI data) should be useful for improving the measurement of 5 and 

^esc,0- 

• Of course, assuming Eqs. ([5]) and (fT0|) we could simply fit for a, S and t csc ,o- However, 
before the role of systematic errors is clarified by a careful analysis, the precise statistical 
significance of such a fit is unclear and we postpone it to a future study. 

For illustration, in Appendix [Al we contrast the data with two specific propagation models, 
a diffusion model and a uniform density model (the latter is a generalization of the usual leaky 
box model). For TZ > 10 GV the results for these models in terms of a, S and i C sc,o are 

a ps 0.7, 5 ps 0, iesc.0 ~ 2.7 x 35 Myr (leaky box), 

a Ps 0.5, S Ps -0.3, t cscfi ps 300 Myr (diffusion). (13) 

3 Lessons from the PAMELA positron data 

Recently, the flux of cosmic ray antiprotons [311 I3"2] and positrons [3J [25] was measured by the 
PAMELA satellite-borne experiment o While the antiproton measurements stand in accord with 
common diffusion models, the positron measurements do not, revealing an excess with respect to 
model predictions. The 'positron excess' triggered numerous attempts to explain the acclaimed 
'non-secondary' origin of the positrons. However, a careful examination of the positron data 
does not reveal any excess with respect to model independent estimates, assuming secondary 
production 22]. Below we analyze the implications of the measured positron fraction under the 
secondary production hypothesis. We also briefly reiterate the main points of [22] . 

A main point we wish to convey is that the PAMELA positron measurement is, first and 
foremost, a measurement of the radiative cooling. To see this, note that it is possible to test 
the secondary production spectrum of the positrons indirectly, by examining the flux of another 
secondary CR which does not suffer propagation energy losses and which originates from the 
same interactions. The first candidate for such a test is the antiproton flux. 

In the left panel of Fig. [3] we show the results of applying Eq. ([1]) to the calculation of the 
antiproton to proton flux ratio. The details of the calculation are given in [22]. We show a 
40% error estimate, representative of relevant cross section uncertainties, as well as the latest 
PAMELA results [32] . Having verified the applicability of Eq. (JTJ for antiprotons, we can now 
find the cooling suppression for positrons by computing the positron flux from Eq. ([1]) , neglecting 
radiative losses, and comparing the result to the observed flux. 

To put the positron and antiproton analyses on equal grounds we give the flux ratio of 
positrons to their leading progenitors, protons. No published PAMELA results currently exist 
for the individual e + ,e~,e ± or p fluxes. Extracting the positron to proton ratio from the re- 
ported e + /e fraction [2 5) requires combining measurements from different experiments. Here 
we perform this combination aiming to illustrate a representation of the data which - unlike 
the e + /e ± fraction - naturally lends itself to theoretical interpretation. A sounder experimental 
procedure will be preferable once complete data are published by the same experiment. We 
interpolate the total electronic flux from FERMI 6 a , available down to TZ = 7 GV, using solar 
modulation parameter $ = 500 MeV. The proton flux we use is that of [M] , for which we adopt 
a 10% normalization error. 

The result is shown in the right panel of Fig. [3] The error band around the theoretical 
upper bound is explained below. As indicated by the arrows, PAMELA measured the cooling 
suppression of the positrons. 

2 See |331 1251 for cautionary notes regarding systematic uncertainties in the positron measurement. 
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Next, we quantify the effect of losses. The positron cooling suppression factor can be defined 
in direct analogy to the suppression due to decay in radioactive nuclei. We find 

f a e+ = ^ « 0.6 x 103 (*\ 05 x J -^, (14) 
' ^Q e+ ^c S c \10GVJ J p (ft)' 1 ' 

for the cross section parameterization of [35] . Using the cross sections from [33] yields a factor 
of ~ 2 increase in the value of f se +, compared to (|14p . for 1Z < 100 GV [37]. This factor of 
two ambiguity is responsible for the main uncertainty in the derivation of the upper bound. 
Considering this ambiguity we adopt a factor of two error estimate, shown in Fig. [3] for the 
theoretical calculation of the positron to proton ratio. Other uncertainties of order 20% arise 
e.g. from the effect of heavy nuclei in the primary CR flux and in the ISM [38 . 

The values of / Sje + for positrons and f a i for some radioactive nucleus should be comparable, 
if examined at the particular rigidity, wherein the corresponding cooling time and observer frame 
lifetime are similar. This observation applies if the cooling time can be approximated by a power 
law, 

t c <xK- Sc , (15) 

and is explained in |22j by noting that decay of radioactive nuclei and energy losses of positrons 
are represented by formally similar terms in a general transport equation. Thus, we expect 

fs, e+ *(-r-)\ (i6) 

\ '•CSC / 

where a and t csc are the same as in Eq. ([5]). Eq. (fTo) is accurate up to corrections of order 
(7 — 1 — S c )~ a , where 7 ~ 3 is the observed spectral index of the positrons. Using a ~ 0.7, in 
accordance with Sec.[H these corrections are at the level of 10% for 5 C ~ 1 (40% for 5 C ~ 0). 
The positron cooling time can be estimated by 

where Ut is the total electromagnetic energy density in the propagation region, averaged over 
CR histories. Comparing f St i of some nucleus with f Syt ,+ of positrons at the same rigidity TV we 
expect 



f S ,e+(K') 



TV \ ( U T 



1.5Myr/ V20GV/ VleVcm 



(18) 



If Ut is rigidity independent, then the ratio Eq. (|18p should rise as 

fs,i/f s ,e+ ~ K 2a . Note that 

a rigidity dependent Ut is possible. In particular, if Ut decreases with increasing CR rigidity, 
then the ratio Eq. (JlSj) could have a milder rise, f s ,i/f s ,e+ °c Tl a ( 1+S ^ with 5 C < 10 

In the rest of this section, we use Eqs. (|16M18[) to analyze the implications of having radioactive 
nuclei data in conjunction with positron data in the same range in rigidity. First, we show that the 
measured positron flux is consistent quantitatively with secondary production. Then, we assume 
a secondary origin for the positrons and derive constraints on the mean traversed electromagnetic 
energy density, Ut- Finally, we comment on the energy dependence of the cooling suppression 
factor of the positrons. 



3 For ~ 50 GV positrons, Klein-Nishina corrections imply that the starlight component in Ut induces <5 C < 1. 
Without extreme assumptions regarding the starlight intensity 1391 . these corrections are of the scale of current 
uncertainties and we defer their discussion to a more refined treatment. 
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Using Eq. (fill) we plot f a e + vs. rigidity in Fig. |H together with the suppression factors due 
to decay for 10 Be, 26 Al and 36 C1 nuclei. Nuclei lifetimes in Myr are specified next to each data 
point. The shaded band denotes the theoretical uncertainty. Error bars show error estimates 
from the experimental procedure. 

Fig. 2] is in broad agreement with the observation made in [22] > namely that the cooling 
suppression factor measured by PAMELA for positrons is consistent quantitatively with that 
found for radioactive nuclei. This observation makes stronger the case for a secondary origin for 
the positrons. In particular, in [22] it was pointed out that f a<e + ~ f s , 10 Be at TZ ~ 20 GV, in 
agreement with Eq. (|18[) for plausible Ut ~ 1 eV/cm 3 . We comment that while the central data 
do not show a clear rising pattern of f s ,i/ f s ,e+ with rigidity, f s .i/ f s .e+ ^ 7?. or so (with, e.g., 
a ~ 0.5 and constant Ut) is perfectly possible within the calculation uncertainties and given the 
visible scatter. 

Assuming secondary positrons, we can use Eq. (|18p together with Fig. 2] to infer the value of 
Ut- First, assuming that the PAMELA result is not grossly wrong, we can set f s e + = f s s&c\ to 
derive an upper bound 



valid in the range TZ ~ 10-40 GV. The bound Eq. ([19]) is nontrivial: for TZ = 40 GV it implies 
Ut < 1.25 eV/cm 3 . For reference, the local value of Ut is determined from CMB, IR and optical 
surveys [ID] to be Iodise = 0(1) eV/cm 3 . Second, Be and Al data give potentially stronger, 
though less robust, constraints 



valid in the same range as Eq. (IT91 . 

Note that Ut < Ucmb ~ 0.25 eV/cm 3 cannot occur physically. If, using Eq. (fT5)l with 
upcoming experimental data, such value will be reliably deduced for Ut it will falsify the as- 
sumption of a secondary origin for the positrons. From Fig. @] we see that CI data at ~ 50 
GeV/nuc will be very useful for such a test. 

Finally we comment on the implications of flat or mildly rising f se + , as derived in Fig. [3] 
from the various data [23 El [33] . By Eq. ([ID]) , the positron suppression factor should rise if t esc 
drops faster than t c with increasing rigidity [22) . The only information which directly constrains 
t csc cx TZ S is the radioactive nuclei data analyzed in Section [2[ This analysis shows that values of 
5 < — 1 are in some tension with observations. However, improved charge ratio measurements, in 
particular Cl/Ar data extending to TZ ~ 100 GV (and so to lifetime tci ~ 40 Myr, comparable to 
the currently existing Be/B measurements), will be required to conclusively rule out (or measure) 
S < — 1 . Regarding the cooling time t c cx TZ~ Sc , effective 5 C < 1 arise if Klein-Nishina corrections 
to Eq. (|17p are important, or if Ut decreases with increasing rigidity. A rigidity dependent Ut 
is observationally and theoretically conceivable. In fact, such behavior of Ut will resemble the 
decreasing traversed matter density inferred from stable secondary to primary nuclei flux ratios. 
To conclude, at present we cannot yet rule out the possibility that 5 C + 5 < 0, for which f se + 
rising with rigidity occurs with secondary production. In the next section we summarize methods 
by which the secondary origin hypothesis for the positrons can be challenged. 

4 Summary and discussion 

In this paper we presented a model independent analysis of the propagation time scales of cosmic 
ray (CR) nuclei and positrons. Motivated by the prospects for improved measurements in the 




(e+ and CI), 



(19) 




(c + , BeandAl), 



(20) 
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Figure 4: The suppression factor due to cooling (e + ) and radioactive decay (nuclei). Numbers 
near nuclei data points denote the lifetime at that point in Myr. 
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near future [H [5J |H1 |S], our analysis demonstrated that it is possible to extract significant 
quantitative information from CR measurements under general assumptions, without committing 
to any particular propagation model. 

In Sec. [2] we studied the radioactive nuclei charge ratio measurements, first analyzed in (11) . 
Examining the suppression of the flux due to decay vs. observer frame lifetime, the clear correla- 
tion across nuclei specie, seen in Fig. (TJ is indicative that the charge ratio analysis is not plagued 
with large systematics. (See, in contrast, [151 [TB].) In addition, the effects of a possible primary 
component in the 26 Al or 36 C1 flux appear to be small on the scale of the measurement errors. 
This makes a combined analysis of the charge ratios a useful model independent tool, first in 
measuring the effect of decay and, second, in constraining the CR residence time. Using charge 
ratio data at 72. ~ 10-40 GV, we showed that a rigidity independent residence time is consis- 
tent with observations, and that the suppression of the flux for all nuclei species is adequately 
described by f s .i — (ti/t eac ) a , with t csc = 100 Myr and a = 0.7. 

In Sec. [3] we analyzed the PAMELA positron data [HH5J. We argued that the positron flux 
is consistent with a secondary origin (22) . Assuming secondary positrons, we derived an upper 
bound on the mean electromagnetic energy density traversed by the positrons, Ut < 1-25 eV/cm 3 
at a rigidity of 1Z = 40 GV. This bound is close to the locally measured electromagnetic energy 
density [3D]- A stronger, though less robust bound can be inferred from Be and Al measurements, 
and may indicate that the positrons spend significant fraction of their confinement time in regions 
far from the Galactic disc, where the electromagnetic energy density is lower than the locally 
observed value. 

Methods by which the secondary origin hypothesis for the positrons can be challenged are 
summarized as follows: 

1. A simple first test for secondary positrons is the condition 



This test is independent of radioactive nuclei data. As seen in Fig. [21 this condition is 
satisfied by current data (22) . 

2. A second, and potentially stronger test for the origin of the positrons is the requirement 



This test can challenge secondary production when f s e+ < 1, and can be performed ro- 
bustly with CI measurements. At present, our evaluation shows that Ut ~ leV/cm 3 > 
UcmBi and so the test is passed in the range where e + and CI data coexist. 

3. Finally, improved radioactive nuclei data will yield a model independent measurement of 
a, describing the functional form of the suppression of the flux due to decay or loss (see 
Eqs. ((5)) . (fT6)) ) and of the rigidity dependence of the residence time, t esc oc 1Z S . Given such 
a measurement, the detailed energy dependence of the ratio f s ,i/f S e+ °£ TZ a ( 1+Sc \ for any 
radioactive nucleus i will further determine the energy dependence of the cooling time, 
t c oc TZ~ Sc . If future measurements (and future PAMELA re-analyses [35]) will confirm a 
rising positron fraction, then a violation of the condition 



fs.e+ < 1. 



(21) 



Ut > Ucmb- 



(22) 



S + 5 C < 



(23) 



will establish a non-secondary origin for the positrons. 
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The effect of decay of radioactive nuclei, defined in Eq. ([5]) and specified by a measurement 
of a, has a deep connection with propagation. In the general class of models in which CR 
production occurs mainly in the thin Galactic disc, the quantity f s> i in the limit of strong losses 
is given by [I] 

fs,i ~ j^f- (24) 

£esc L'i 

Here U is the observer frame lifetime, Li is the scale height above the disc occupied by the 
decaying nuclei and t esc is the escape time from the propagation region of scale height L csc . 
t csc and L csc are defined for stable CRs, and would apply to the specie i had its lifetime been 
infinite [13] ■ The functional form of f s< i is determined by the nature of the CR motion in the 
trapping volume. Diffusion models have L ~ \ft and give f S; i ~ \/ti/t csc or a ~ 0.5. Leaky 
box models impose L csc — Li, giving / S) j ~ ti/t esc or a ~ l|j Models of compound diffusion [T] 
predict L ~ t 25 and thus f Sj { ~ (ti/t csc ) - 75 or a ~ 0.75. It is worth pointing out in this 
context that the currently fashionable homogeneous diffusion models with large halo and thin 
disc predict a = 0.5 to high accuracy, see Appendix [A] This value of a is currently allowed by 
data. Establishing a > 0.5 will rule out this class of models. 

The AMS-02 experiment [5] has the capabilities to take the analysis to a new level of accuracy 
in the near futureO AMS-02 will resolve CR charge all the way to Z — 26 and up to TeV/nuc. In 
addition, AMS-02 will measure isotopic ratios including 9 Be/ 10 Be up to 10 GeV/nuc [32]. This 
measurement will be a valuable check of the interpretation of the charge ratios. Finally, AMS-02 
will also measure the antiproton and positron fluxes up to ^300 GeV [5] , assessing the validity 
of the PAMELA measurements [3T] [32] HI ES] and putting the combined analysis of positrons 
and radioactive nuclei on firmer observational grounds. 

Further progress which is potentially within immediate reach will follow a publication of the 
individual spectra of e + ,e~,p and nuclei by the PAMELA collaboration. These data can be 
used to extract quantities such as e + /p and e + /p, which are theoretically clean probes of CR 
propagation, in a direct manner. 
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A Propagation model examples 

Here we demonstrate the implications of the radioactive nuclei data for two specific models. First 
we consider a version of the leaky box model, with spatially uniform CR distribution at a given 
rigidity but with a rigidity dependent halo scale height L = L(1Z). Second we consider thin disc 
diffusion models with a large halo. The models predict: 

fs,i = , 1 , , leaky box, 

-L T" £csc / H 

fs,i = \At/*esc tanh ( v^cscAi) , diffusion. (25) 



4 As shown in AppcndixfX] for the leaky box model the strong loss result ceases to apply for ti/t eBC > 0(0. 1), 
leading to a ~ 0.7 < 1. The strong loss result a ~ 1 could in principle be tested using short-lived nuclei such as 
I4r< 

6 °- 

5 We note that in |41| the implications of future AMS-02 measurements were analyzed in the context of a 
diffusion model. 
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We fit the suppression factor for 10 Be, 26 Al and 36 C1 using Eqs. (j2"5")) and (fTU|) . so that our fit 
parameters are i osc ,o and S. In each case we arrange the parameters of the model (halo size and 
escape time for the leaky box, and halo size and diffusion coefficient for the diffusion model) 
such as to reproduce the measured CR grammage 29J. The results for both models are shown 
in Fig. [SJ We make the following comments: 

• The suppression factor / S) j in Eq. (|25p involves the nucleus specie label i only via the 
lifetime ij . Any explicit dependence on the fragmentation cross sections and CR grammage 
is canceled. This would not have been the case had we worked with the surviving fraction 
fi (as easily verified from Eq. ©). We expect this result to hold in general for disc+halo 
models, in which the CR composition in the regions where the secondaries are produced is 
similar to the CR composition near earth. 

• For 1Z > 10 GV, values of f s .WBc, fs, 26 Ai and f s , 3e ci obtained from Eq. (|25p satisfy the 
power law assumption Eq. to better than 10%, with a =0.7 (0.5) for the leaky box 
(diffusion) cases. For the leaky box model, the parameter t esc has the meaning of 2.7 times 
the escape time of the CRs. For the diffusion model, t csc is the average time it takes a CR 
produced on the disc to reach the halo boundary. 

• Both the leaky box and diffusion models accommodate the data reasonably well. We find 

t QSC « 2.7 x (30- 45) x (ft/10GV) 0±0 ' 2 Myr, leaky box, 
t csc « (200 - 400 ) x (ft/10 GVr°- 3±0 - 3 Myr, diffusion. 
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